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Preface 

The  use  of  a one-dimensional  model  atmosphere  calculation  incorporating 
dissociation  of  O.,  and  ionized  species  demonstrates  that  dissociation,  through  the 
production  of  atomic  oxygen,  can  drive  the  vertical  thermosphere  and  upper 
mesosphere  into  diurnal  oscillations.  It  is  further  demonstrated  that  neutral 
mesospheric  and  lower  thermospheric  turbulence,  through  its  control  of  atomic 
oxygen  flow,  can  cause  large  amplitude  variations  of  these  oscillations. 
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Dissociation 


Driven  Diurnal  Oscillations 


I.  I INTRODUCTION 


The  dynamic  processes  taking  place  among  the  neutral  species  in  the  thermos- 
phere are  of  an  extremely  complex  nature.  The  interplay  of  driving  forces  and 
atmospheric  response  has  been  considered  in  many  theories  and  texts  (to  cite  a 
few)  ’ ' ' and  many  more  in  the  recent  literature.  Generally,  heat  energy  is  con- 
sidered as  the  driving  mechanism  of  tides  and  motions  in  these  theories.  In  the 
thermosphere,  the  EUV  solar  radiation  below  1000  A is  absorbed  above  120  km  and 
heats  the  neutral  gas  through  ionization,  dissociation,  and  recombination  of  the 
plasma.  Radiation  at  longer  wavelengths  is  absorbed  below  120  km  and  again  there 

is  the  generation  of  heat  energy  through  absorption,  dissociation,  and  subsequent 

4 5 

recombination  of  ground  state  and  excited  species  ’ that  drives  the  atmosphere  to 


(Received  for  publication  30  May  1978) 

1.  Chapman,  S.  , and  Cowling,  T.G.  (1952)  The  Mathematical  Theory  of  Non- 

Uniform  Gases,  2nd  ed.  Cambridge  University  Press,  New  York. 
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3.  Harris,  1.,  and  Priester,  W.  (1962)  Time-dependent  structure  of  the  upper 

atmosphere,  J.  Atmos,  bci.  19:286-301. 
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4.  Lindzen,  R.  W. , and  Ulake,  D.  (1970)  Internal  Gravity  waves  in  atmosphere  with 

realistic  dissipation  and  temperature,  2,  Geophvs.  fluid.  IJyn.  ,2:31 

5.  Volland,  M.  , and  Mayr,  H.  G.  (1973)  A numerical  study  of  three-dimensional 

variations  in  the  thermosphere,  Ann.  Geophys.  29:1  102. 
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respond  to  this  diurnal  stimulus.  We  would  like  to  discuss  another  mechanism 
that  can  induce  atmospheric  motions.  It  is  ignored  in  most  calculations  but  when 
incorporated  it  is  generally  not  considered  separately  from  the  heat  source.  It  is  usually 
taken  into  account  in  one-dimensional  models  that  try  to  display  minor  species  and 
temperature  distributions  and  diffusive  motions  in  the  lower  thermosphere  and 
ionosphere.  This  mechanism  is  the  dissociation  of  molecular  oxygen  to  produce 
oxygen  atoms  that  through  their  creation  and  subsequent  diffusive  migration  from 
their  source  will  force  the  atmosphere  into  motion.  This  effect,  which  in  our  model 
exhibits  a diurnal  response  to  the  diurnal  driving  force,  is  created  when  there  is 
complete  coupling  among  all  species  through  collisions.  Interactions,  primarily 
atomic  oxygen  with  Ngr  start  the  molecular  nitrogen  into  a diurnal  up  and  down 
motion.  This  action  continues  until  all  species,  neutral  and  ionic,  react  to  this 
mean  mass  motion  through  collisions  or  through  chemical  interactions. 

The  theory  we  apply  follows  from  Chapman  and  Cowling,  (Reference  1)  and 

g 

Hirschfelder  et  al  for  the  equation  of  motion  in  a coupled  multicomponent  gas 

7 

(also  see  Banks  and  Kockarts  ).  Expressed  as  one-dimensional  equations  in  the 
vertical,  thera^ire  the  neutral  motions 
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In  these  equations 


C.  is  the  velocity  of  the  i th  species, 
n.  is  the  concentration  of  the  i th  species. 


n 


is  the  electron  concentration, 

T is  temperature, 

k is  Boltzmann's  constant. 


6.  Hirschfelder,  J.O.,  Curtiss,  C.  F. , and  Bird,  K.  B.  (1954)  Molecular  Theory 

of  Gases  and  Liquids,  John  Wiley,  New  York. 

7.  Banks,  P.  M.  , and  Kockarts,  G.  (1973)  Aoronomy,  Academic  Press,  New  York. 
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is  the  collision  frequency  of  species  i in  species  j,  = Zy 

is  the  gyrofrequency  of  the  electrons, 
i.3  the  magnetic  dip  angle. 


These  equations  are  coupled  to  the  equation  of  continuity 


9 n.  9 0. 

W1  = Fi  ' ni  Ri  ' TT 

Dp.  C. 
i i i 

and  the  mean  mass  motion  (C  ) is  determined  from  C = 

o O V , 

DPi 

i 1 


(3) 


where 

F . is  the  rate  at  which  epecies  i is  chemically  formed, 

n.K.  is  the  rate  at  which  species  i is  chemically  removed, 

0.  = n.  CL  is  the  flux  of  species  i, 

p^  is  the  mass  density  of  species  i. 

The  initial  conditions  used  in  these  calculations  are  the  results  discussed  in 
Philbrick  et  al  that  displays  a one-dimensional  model  of  the  atmosphere  using  the 
temperature  and  mass  density  data  as  well  as  the  turbulent  diffusion  coefficients 

8 

(K  and  3 K)  determined  from  the  ALADDIN  1 experiments.  These  Philbrick  et  al 
data  were  then  extrapolated  from  150  km  to  400  km  using  the  diffusion  equilibrium 
hypothesis.  Further  results  giving  the  minutiae  of  the  expanded  ion-neutral  model 
synopsized  here  will  be  presented  in  a more  complete  dissertation  later.  For  this 
report,  we  shall  present  only  some  of  the  salient  points  supporting  our  thesis  from 
the  fifteenth  model  day  when  most  of  the  major  species  have  arrived  at  near  diurnal 
reproducibility.  These  calculations  cover  the  altitude  region  50  to  400  km. 


2.  RKSUI.TS  AND  CONCLUSIONS 


The  results  of  the  calculation  that  display  the  diurnal  forcing  by  the  oxygen  flow 
that  in  turn  drives  the  other  constituents  are  presented  in  Figures  1,  2,  3,  and  4. 
These  show  the  diurnal  variation  of  O,  Ar,  N^.  and  p at  specific  altitudes  in  the 
mesosphere  and  thermosphere  for  the  turbulent  diffusion  profile  (K)  given  in  Fig- 
ure 5.  It  is  obvious  that  the  diurnal  variations  of  the  neutral  constituents  peak 

8.  Philbrick,  C.  K. , Narcisi,  11.  S.  , Good,  ILL.,  Hoffman,  li.  S.  , Keneshea,  T.J. 
MacLeod,  M.A.,  Zimmerman,  S.  P.  , and  Heinisch,  B.  W.  (1973)  The 
ALADDIN  experiment-Part  II,  Composition,  Space  lies.  13:441-448. 
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at  ~ 1400  hr,  a phase  time  in  accord  with  the  recently  reported  data  of  ESHO-4.  ' 
Indeed,  the  comparison  for  the  O/N^  ratio  at  250  km  (f  igure  6)  demonstrates  fair 
agreement  in  amplitude  and  excellent  agreement  in  phase  between  theory  and  experi- 
ment. In  tlie  mesosphere,  on  the  other  hand,  there  is  a reversal  in  the  phase  of  O 
as  compared  with  that  in  the  thermosphere,  particularly  observed  around  the  95-km 
peak.  Here  the  maximum  occurs  at  ~0200  hr,  in  phase  with  the  maximum  of  the 
northern  hemisphere  midlatitude  green  line  measurements  (Figure  7)^  where  the 
relative  model  variation  in  O is  ~ 15  percent,  similar  to  that  determined  from  green 
line  (0,  5577  A)  variations. 


LOCAL  TIME  (hr) 


Figure  1.  The  Diurnal  Variation  of  Atomic  Oxygen 
Concentration  for  Selected  Altitudes  C ling  the  Tur- 
bulent Diffusion  Coefficients  Labeled  K.  Note  the 
reverse  behavior  of  the  concentration  at  95  km 
compared  to  that  of  105  km 


9.  Kohnlein,  W.  , Trinks,  IF,  and  Volland,  H.  ( 1975)  The  O and  Ng  density  ratio 
in  the  thermosphere  derived  from  FSRO-4  data.  Space  Kes.  15: 207 -291. 

10.  Brenton,  J.  G.  , and  Silverman,  S.  M.  (1970)  A study  of  the  diurnal  variation  of 
the  5577  A OI  airglow  emission  at  selected  1GY  stations,  Planet.  Space  Sci. 
18:641-653.  
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Figure  6.  The  O/N2  Ratio  at 
Thermospheric  Altitudes  and 
Compared  to  the  Measurements 
of  Kohnlein  et  al® 
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Figure  7.  The  Nighttime  Variation  of  the  Integral 
of  the  Atomic  Oxygen  Concentration  (dashed  line) 
Where  the  Limits  of  Integration  Are  F rom  85  to 
120  km.  The  ground  based  green  line  (5577  X) 
measurements  (solid  line)  at  Sac  peak  (Brenton 
and  Silverman)1®  are  compared  to  O column 
density  ( / O^  dz)  with  good  agreement 


The  mean  mass  motion,  defined  previously,  follows  the  above  picture,  upward 
during  the  dissociation  interval,  maximizing  at  the  time  of  the  peak  of  the  thermos- 
pheric density  variations,  reversing  after  sunset  and  achieving  maximum  downward 
amplitude  at  the  time  the  altitude  invariant,  atomic  oxygen  peak  at  95  km  maximizes 
These  results  combine  to  form  the  picture  initially  stated  wherein  the  solar 
dissociation  flux  provides  the  driving  mechanism  to  set  the  atmosphere  into  diurnal 
oscillation  without  additional  forces.  The  amplitude  of  this  pressure  force,  deter- 
mined in  this  calculation  due  to  dissociation,  is  approximately  10  percent  of  the 
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solar  thermal  driving  force,  peaking  at  ~ 130  km.  The  above  action  being  achieved 
primarily  through  the  interaction  of  the  oxygen  atoms  with  all  other  atmospheric 
species.  If  this  picture  is  correct,  then  a variation  of  the  turbopause  or  turbulent 
intensities  should  have  a significant  effect  upon  the  diurnal  oscillation  of  the  atmos- 
phere. 

This  we  shall  show  is  the  case  following  Keneshea  and  Zimmerman1'  who  have 
demonstrated  that  the  flux  of  atomic  oxygen  exhibits  a strong  diurnal  variability 
when  the  turbopause  is  near  or  lower  than  95  km.  However,  when  the  turbopause 
was  raised  to  112  km,  their  results  demonstrated  that  the  diurnal  variability 
disappears,  and  only  downward  motion  of  G prevails.  If  the  above  is  correct,  then 
the  driving  mechanism  of  the  diurnal  oscillations  hypothesized  here  will  be  minimized 
or  disappear  completely  when  we  increase  the  turbulent  diffusion  coefficient  and/or 
raise  the  turbopause.  With  this  in  mind,  we  now  examine  the  results  for  calcula- 
tions using  turbulent  diffusion  coefficients  a factor  of  3 larger  than  those  shown  in 
Figure  5.  In  essence,  using  3K  has  the  effect  of  both  raising  the  turbopause  and 
increasing  the  turbulent  intensity.  The  data  examined,  but  not  displayed  here,  is 
that  of  modal  day  21,  and  although  diurnal  reproducibility  has  not  yet  been  achieved, 
nonetheless,  the  results  using  3K  are  extremely  different  from  those  using  K.  The 
diurnal  oscillations  using  3K  are  much  smaller,  indeed  practically  non-existent  and 
while  the  phase  of  the  thermospheric  oscillat  ms  are  the  same  as  those  uaing  K, 
the  mesospheric  variations  of  O are  now  reversed.  This  puts  the  mesospheric  and 
thermospheric  density  variations  in  phase  so  that  the  atomic  oxygen  peak  now  maxi- 
mizes at  ~ 1400  hours.  This  then  would  cause  the  nighttime  variation  of  green  line 
to  exhibit  a minimum.  Figure  8 shows  that  indeed  there  is  data,  measured  in  the 
southern  hemisphere111  in  which  the  nighttime  intensity  of  green  line  does  minimize. 
The  mean  mass  velocity,  while  still  showing  a diurnal  oscillation,  is  significantly 
smaller  in  amplitude  than  that  using  K.  These  results,  coupled  to  those  of  Keneshea 
and  Zimmerman11  confirm  the  hypothesis  set  forth  here  that  the  diurnal  variation 
of  the  atomic  oxygen  flux,  and  its  control  by  turbulence  distribution  and  intensity 
in  the  lower  thermosphere  and  upper  mesosphere  indeed  does  add  to  the  control  of 
the  mass  motion  in  the  thermosphere  and  into  the  mesosphere. 


11.  Keneshea,  T.  J.  , and  Zimmerman,  S.  P.  (1970)  The  effect  of  mixing  upon 
atomic  and  molecular  oxygen  in  the  70-170  km  region  of  the  atmosphere, 
J.  Atmos.  Sci.  27:831-840. 
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Figure  8.  The  Nighttime  Variation  of  Green  Line  Measurements 
at  Various  Selected  Sites  Demonstrating  the  Nighttime  Minimum 
Corresponding  to  a Minimum  in  / dz 


3.  CONCLUSIONS 

We  have  shown  that  a somewhat  significant  driving  force  of  thermospheric  and 
upper  mesospheric  diurnal  variations  is  the  production  of  atomic  oxygen  by  photo- 
dissociation and  its  redistribution  under  the  control  of  local  gradients  due  to  the 
interaction  of  production,  molecular  diffusion,  and  turbulence.  For  the  turbulent 
condition  that  creates  the  turbulent  diffusion  coefficient  K,  we  observe  that  varia- 
tions of  O,  Ng,  Ar,  and  follow  well  the  measurements  of  these  species  in  the 
thermosphere  obtained  from  satellite  mass  spectrometer.  Further,  the  nighttime 
intensity  variations  of  the  green  line  can  be  approximately  reproduced  in  amplitude 
and  phase  following  the  O variations.  When  the  turbulent  motions  are  intensified, 
the  atomic  oxygen  flux  is  directed  more  into  the  downward  direction,  and  thus  the 
upward  forcing  of  the  diurnal  oscillations  is  significantly  reduced.  Also,  green 
line  variability  is  slightly  dampened  and  reversed  in  phase  from  that  due  to  the 
smaller  turbulent  diffusivity.  Further,  we  see  that  not  only  dees  turbulence  control 

12.  Zimmerman,  S.  P.  , and  Trowbridge,  C.  A.  ( 1973)  The  measurement  of 

turbulent  spectra  and  diffusion  coefficients  in  the  altitude  region  95  to  110  km. 
Space  Res.  13:203-208. 
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the-  atomic  oxygen  flow  but,  through  interactions  in  the  form  of  collisions,  other 
thermospheric  species  are  induced  into  synchronous  motion  until  the  total  mass 
density  is  forced  into  diurnal  oscillations.  However,  a more  accurate  picture  will 
be  derived  with  the  inclusion  of  thermal  forcing  in  a three-dimensional  calculation. 
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